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CN| , The production of different particle species is recently measured in Pb — Pb collisions by the 

ALICE experiment at ^/s = 7 TeV. This motivates the use of various bosons and baryons measured 
■ at lower center-of-mass energies in comparing their ratios to the hadron resonance (HRG) gas model 

04 \ and PYTHIA event generator. It is found that the particle-to-antiparticle ratios are perfectly 

^ ■ reproduce by means of HRG and PYTHIA at RHIC and LHC energies. The kaon-to-pion and 

^ ■ proton-to-pion ratios are entirely overestimated by the HRG model. The PYTHIA event generator 

I obviously underestimates the kaon-to-pion ratio and simultaneously reproduces the proton-to-pion 

ratio, almost perfectly, especially at LHC energy. While matter-to-antimatter and non-strange 
abundances are partly in line with predictions from the HRG model, it is found in the ALICE 
experiment that the measured baryon ratios are suppressed by a factor of ~ 1.5. The strange 
abundances are overestimated in the HRG model. 
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I. INTRODUCTION 



^ \ Understanding the dynamical properties of hot and dense hadronic matter is one of the main motivations of 

0^ . heavy-ion experiments, whieh in turn offer unique possibilities to study hadronic matter under extreme condi- 
I tions [H-Q and to compare with the lattice QCD simulations Q. The Relativistic Heavy- Ion Collider (RHIC) 
. has shown that the bulk matter created in such collisions can be quantitatively described by hydrodynamic 
models The created hot and dense partonic matter, the quark-gluon plasma (QGP) or colored quarks and 
gluons, where quarks and gluons can move freely over large volumes comparing to the typical size of a hadron, 
rapidly expands and cools down. Over this path, it likely undergoes phase transition(s) back to the hadronic 
matter. Different thermal models can very well reproduce the produced particle abundances, which are governed 
in chemical equilibrium by two parameters, the chemical freeze-out temperature Tch and the baryochemical po- 
tential /Zb, where the latter reflects the net baryon content of the system, directly, and the center-of-mass energy, 
indirectly. 

5^ ■ That the particle abundances at Alternating Gradient Synchrotron (AGS), Super Proton Synchrotron (SPS) 

! and RHIC energies are consistent with equilibrium populations @ makes it possible to extract both freeze- 
out parameters over a wide range of center-of-mass energies ^/snn from fits of measured particle ratios with 
thermal models, which obviously do not count for hadrons interactions in the final state. Nevertheless, the 
formation of resonances can only be materialized through strong interactions since the resonances (fireballs) 
are composed of further resonances (fireballs), which in turn consist of resonances (fireballs) and so on 0]. 
Taking into consideration all kinds of resonance interactions by means of the S-matrix, which describes the 
scattering processes in the thermodynamical system, reduces the resulting virial term, so that the partition 
function turns to be reduced to the non- relativistic limit, especially at narrow width and/or low temperature 
T . The resonance contributions to the partition function are the same as that of collisionless particles with 
some effective mass. All possible interactions modifying the relative abundances are found negligible in the 
hadronic phase d, Q. It is assumed that the hadron resonances with masses < 2.8 GeV avoid the singularities 
expected at Hagedorn temperature [lol - [l3j . 

The earliest idea about enhancement of strangeness abundances as a key signature of QGP formation in 
heavy-ion collisions has been proposed, three decades ago [l^ . During the hadronization process, T drops from 
about two hundred to few tens MeV, the dynamical mass of strange quark drops, as well. As a consequence. 



'Electronic address: |a. tawfik@eng.mti.edu. eg) Electronic address: Iatawfik@cern.ch1 



2 



strangeness in QGP would equilibrate on small time scales relative to those in a liadronic matter [l5[. At SPS 
energy, there were strong enhancements observed in AA- and almost neglecting ones in pA-collisions. This 
has been cited as an experimental evidence for QGP formation [l6j . The enhancements at AGS energy are 
considered a typical rescattcring of produced hadrons. Assuming that thermally equilibrated QGP hadronizes 
into a maximum entropy state, a test for strange quark saturation in the early stages is provided by comparing 
final state hadron yields to thermal model predictions [13, [ill ■ 

Studying the ratios of particle yields helps in determining the freeze-out parameters and in eliminating the 
volume fluctuations. Furthermore, the dependence of the freeze-out surface on the initial conditions can be 
neglected. Using statistical model with corresponding choices of the thermal model parameters with ^^snn 
made according to the systematics extracted from h eavy -ion collisions at lower energies, some predictions of 
particle abundances at LHC energy are reported Using the hadron resonance gas (HRG) model, p/p 

and K~ /K'^ are given in dependence on ^/Tnn j2lll22l|. The present work is motivated by recent measurements 
of pion, kaon and proton ratios in central Pb-Pb collisions at ^/sNN = 2.76 TeV by the ALICE experiment. 

The present paper is organized as follows. Section|n]elaborates details on HRG model and the event generator 
PYTHIA. The results and discussions are outlined in section Hill The conclusions are discussed in section HVl 

II. THE MODEL 

A. Hadron Resonance Gas Model 

Treating the hadron resonances as a free gas is supposed to construct the thermodynamic pressure in the 
hadronic phase i.e., < [lol - [l3j . This is valid for coUisionless and interacting hadron resonances. It has 
been shown that the thermodynamics of strongly interacting s yst em can also be approximated to an ideal 
gas composed of hadron resonances with masses < 2 GeV [^, UM- Therefore, the confined phase of QCD, 
the hadronic phase, would be modelled as a non-interacting gas of hadron resonances. The grand canonical 
partition function can be given as: 

Z{T,V) = Tr exp""/"^] , (1) 

where H is the Hamiltonian of the system, which is given by the summation over the kinetic energies of 
relativistic Fermi and Bose particles. The main motivation of using this Hamiltonian is that it contains all 
relevant degrees of freedom of confined, strongly interacting matter. It includes implicitly the interactions 
that result in the formation of the resonances. In addition, it has been shown that this model gives a quite 
satisfactory description of the particle production in the heavy-ion collisions. With the above assumptions, the 
dynamics of the partition function can be calculated exactly and be expressed as summation over single-particle 
partition functions of all hadrons and their resonances. 

\nZ{T,^,,V) = VlnZi(r,F)=^±^ r edk\n{l±cxp[{^l,^e,)/T]}, (2) 

where £i{k) — (fc^ -f m|)^/2 is the i— th particle dispersion relation, gi is spin-isospin degeneracy factor and ± 
stands for bosons and fermions, respectively. 

Before the discovery of QCD, it was speculated about a possible phase transition of a massless pion gas to a 
new phase of matter. Based on statistical models like Hagedorn [23| and Bootstrap [25|, the thermodynamics of 
such an ideal pion gas is studied, extensively. After the QCD, the new phase of matter is known as QGP. The 
physical picture was that at Tc the additional degrees of freedom carried by QGP are to be released resulting in 
an increase in the thermodynamic quantities. The success of HRG in reproducing lattice QCD results at various 
quark fiavours and masses (below T^) changed this physical picture, drastically. Instead of releasing additional 
degrees of freedom at T > Tc, it is found that the interacting system reduces its effective degrees of freedom at 
T < Tc- In other words, the hadron gas has much more degrees of freedom than QGP. 

At finite temperature T and baryo-chemical potential ^i, the pressure of i-th hadron or resonance reads 

^ a r°° 

p{T,pi,) = ±Y,^T fc2^fcln{l±exp[(/i,-eO/r]}, (3) 
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where N is the total number of hadron resonances of interest. As no phase transition is conjectured in HRG, 
summing over ah hadron resonances results in the final thermodynamic pressure in the hadronic phase. Switching 
between hadron and quark chemistry is given by the correspondence between the hadronic chemical potentials 
and that of the quark constituents, for example, fii ~ Srihfiq + UgfJ-s, where nb{ns) being baryon (strange) 
quantum number. The chemical potential assigned to the degenerate light quarks is fiq = (pu + fJ-d)/'2 and the 
one assigned to strange quark reads /is = /ig — /^s . 

The strangeness chemical potential /ig is treated as a dependent parameter. Basically, it is calculated as 
a function of T and fXb, based on the fact that the overall strange quantum number has to remain conserved 
in heavy- ion collisions Q. Therefore and in order to assure vanishing strange charge, iis{iJLb,T) has to be 
calculated, whenever baryochemical potential ^b and/or temperature T arc changed. 

When ignoring all decay channels, the particle number density is given by the derivative of the partition 
function, Eq. |3j with respect to the chemical potential of interest. 



, , , dkk^ -, (4) 

i 

In the present work, T and [ib, at which the chemical freeze-out takes place, are characterized by constant s/T^ , 
where s is the entropy density. Details about this ratio and its physical meaning are given in [13|. 



B. PYTHIA 



Although the experimental data shown in the present work are taken from heavy-ion collisions, the comparison 
with PYTHIA. which is designed to generate multiparticlc production in collisions between elementary particles, 
e^e~ , pp and ep, is possible at very high energies. The bulk of PYTHIA multiplicities is formed in jets, i.e. in 
collimated bunches of hadrons or resonances decaying into further hadrons produced by the hadronization of 
partons [1^ . The relative proportion of strange particles is as expected smaller in comparison with nonstrange 
hadrons. PYTHIA is capable of simulating for different processes including hard and soft interactions, parton 
distributions, initial/final-state parton showers, multiple interactions, fragmentation and decay. 

The measured particle production is conjectured as an indicator for the formation of QGP, especially in 
heavy-ion collisions. In pp collisions, the spatial and time evolution of the system is too short to assure initial 
conditions required to drive hadronic matter into partonic QGP. The HRG model has been successfully used 
to describe particle abundances and their fluctuations in heavy-ion collisions. The thermodynamics of lattice 
QCD is also reproduced by means of HRG at temperatures below the critical one 

In Refs. [llljll^l, we have noticed that the collective flow of strongly interacting matter in heavy- ion collisions 
makes HRG model underestimating the particle yields ratios measured in pp collisions, especially at low energies. 
Nevertheless, it was found that the differences between particle yields ratios of pp and AA collisions almost 
disappear, at LHC energies. In light of this, the comparison with PYTHIA remains an enlightening feature. 
Although, it gives comparable high-energy results as the ones from heavy-ion collisions, its initial conditions 
would be reflected in the collective properties in the final state. This would include - among other - the issue of 
strangeness suppression, as the mass of strange quark is heavier than that of up and down quarks. Therefore, 
the production of strange hadrons is generally suppressed relative to hadrons containing only up and down 
quarks. 



III. RESULTS AND DISCUSSIONS 



In Fig. [1] the ratios of pion, kaon and protons at RHIC and LHC energies (open symbols) are compared with 
the HRG model (horizonal lines) and the PYTHIA event generator (solid circle). The HRG results are calculated 
at 200 GeV corresponding to the RHIC energy and 2.76 TeV corresponding to the LHC energy. The PYTHIA 
simulations are performed for pp collisions at 2.76 TeV. First, we study the ratios of particle-to-antiparticle 
given in the left panel. It is obvious, that the experimental results are well reproduced by means of both HRG 
and PYTHIA. Furthermore, the energy-dependence is very well reflected. The right panel shows the different 
proton-to-pion ratios. We notice that PYTHIA results describe very well the ALICE experimental results taken 
at 2.76 TeV. It is apparent that the HRG results on p/ir nearly reproduce the experimental data, while the ratios 
oip/ir and (p + p)/{Tr'^ +tt^) obviously underestimate the corresponding ALICE results. The RHIC results on 
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the three ratios are very weh reproduced by means of the HRG modeL It is interesting to notice that the middle 
panel gives amazingly different results. Both experimental data sets are not clearly distinguishable, although 
the huge jump in center-of-mass energy from RHIC (200 GeV) to LHC (2760 GeV). We notice that the 

PYTHIA event generator entirely underestimates the experimental results, while the HRG model overestimates 
them. 
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Fig. 1: Measured particle ratios (symbols) are compared to the predictions from HRG model (horizontal lines) and 
PYTHIA even generator (solid circles). 



Fig. [T] shows a collective comparison. Another comparative study is summarized in Tab. HI where the HRG 
results on particle-antiparticle ratios are compared with the corresponding results reported in Ref. [l9| . The 
mixed particle ratios are listed in Tab. HH It compares the HRG results with previous studies [T9l.[20l|. 





A" 




p 


A 




n+ 




A'+ 




p 


A 




Q- 


HRG 


0.994 


0.9997 


0.978 


0.983 






Ref. [19] 


0.9998 


0.9998 


0.989 


0.992 


0.994 


0.996 



Tab. I: Ratios of matter-antimatter compared with the results reported in Ref. [l^. 
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Tab. 11: Ratios of mixed particle species compared with the results reported in Ref. [T^ and and compared with 
experimental measurements at RHIC and LHC. 

A quantitative comparison is illustrated in Fig. [2l [3] and ID In Fig. [2l tt^/tt^, K^/tt^ and P/ty^ ratios 
in dependence on ^/snn- The solid curve represents the HRG calculations for tt^/tt^. As shown in Fig. [l] 
the agreement with both experimental data sets (RHIC and LHC) is excellent. This has been introduce in 
Ref. [lllll^. The dashed curve represents the energy evolution of K~ j-n^ ratio. Here, the experimental data 
lay below the experimental ones. The third ratio is given by the dotted curve. We find that the HRG model 
reproduces the RHIC results for p/7r+, for which the ALICE results are overestimated. 

In Fig. 131 the strange bosonic ratio of K~ /K'^ , strange to non-strange bosonic ratio of j-n^ ^ and baryon 
to non-strange boson p/7r+ ratios are given in dependence on y^JJvw- The solid curve represents the HRG 
calculations for anti-kaon to kaon ratio K~ jK^ . As shown in Fig. [1] the agreement with both experimental 
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Fig. 2: Ratios of tt'^/tv , K /tv^ and p/yr^ are given in dependence on y^sjvjv and compared with experimental 
measurements at RHIC and LHC. 



data sets (RHIC and LHC) is excellent [2l|, (2^. The particle ratio of A'+/7r+ is represented by the dashed 
curve. It slightly decreases with increasing y/sJjN- Once again, the HRG model obviously overestimates the 
experimental results, entirely. The dotted curve represents the HRG results on p/n^. We find that the RHIC 
results are well described by the HRG model, while the ALICE data are obviously overestimated. 




Fig. 3: Ratios of K /K'^, K'^ /% , and p/vr"*" are given in dependence on ,/snn and compared with experimental 
measurements at RHIC and LHC. 

Finally, the three remaining ratios p/p, {K'^ + K~)/{tt~^ + tt~) and {p +p)/(7r+ + tt~) are given in Fig. H] 
as function of y/SNN- Once again, the agreement between the experimental results (RHIC and LHC) and the 
HRG model on p/p is excellent, Fig. [T]and Ref. [2ll. [22j. The kaon-to-pion ratio is overestimated by the HRG 
model, while the proton-to-pion ratio is well reproduced at the RHIC energy but once again overestimated at 
the LHC energy. The kaon-to-pion ratios are entirely overestimated. Along the entire energy range, the HRG 
model gives almost 30% higher values than the measured ratios. 
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Fig. 4: Ratios of p/p, {K'^ + K~)/ (tt"*" + 7r~), and {p + p)/ (vr"*" + 7r~) are given in dependence on ^/snn and compared 
with experimental measurements at RHIC and LHC. 

IV. CONCLUSIONS 

Based on the study given in section Hill two remarks are now in order. The first one is that the energy scan 
of PYTHIA is completely excluded, because of its basic assumptions. The PYTHIA simulations are performed 
for pp collisions. Therefore, the comparison with the heavy-ion collisions would not helpful, especially at low 
energies. The second remark is that the HRG model counts for all baryonic and bosonic resonances with masses 
< 2.8 GeV, assures conserved strange degrees of freedom, excludes all decay channels and finally characterizes 
the chemical freeze-out by constant s/T^ over the whole energy range. 

The number densities of different particle species are calculated in framework of the HRG model and compared 
with the experimental results at RHIC and LHC energies. Fig. [T]summarizes the results, collectively. We notice 
that the particlc-to-antiparticle ratios are best described by the HRG model and PYTHIA event generator 
at both RHIC and LHC energies. The results on mixed particle ratios are different. Obviously, the strange 
particles are not reproduced by the HRG model. At RHIC and LHC energies all ratios including kaon-mesons are 
overestimated. The PYTHIA event generator apparently underestimates all these ratios, especially at 2.76 TeV. 
The proton-to-pion ratios are partially described by the HRG model, at RHIC energy, while the ALICE results 
are apparently overestimated. Obviously, PYTHIA reproduces very well all these ratios, especially at 2.76 TeV. 

So far, we conclude that the particle-antiparticle ratios are very well reproduced by means of HRG and 
PYTHIA. These two models are partially able to reproduce the mixed particles ratios. This appears very clearly 
when comparing their results with the experimental kaon-to-pion and proton-to-pion ratios. The strange ratios 
quark flavours seem to play an essential role in explaining the discrepancy with the kaon-to-pion ratios. On one 
hand, the HRG model, which excludes all decay channels, seems to overestimate the experimental results. On 
the other hand, the PYTHIA event generator, which is basically originated from pp interactions, whose spatial 
and temporal evolutions are likely short to assure initial conditions driving the hadronic matter to the partonic 
QGP, seems to underestimate the experimental results. The threshold for strange quark production in QGP is 
much smaller than in hadronic matter. This effect would be further enhanced in baryons with multiple strange 
quarks. While QGP is more likely to be found in heavy-ion collisions, the strangeness enhancement in high 
energetic pp collisions would be a sign of a collective effect [l^l . 

So far, there are three main features of ALICE measurements. Fig. [T] The first one is that the particle-to- 
antiparticle we very well reproduced by means of the HRG model and the PYTHIA event generator at RHIC 
and LHC energies. The second one is that the kaon-to-pion and proton-to-pion ratios are entirely overestimated 
by the HRG model. The third one is that the PYTHIA event generator obviously underestimates the first ratio 
and simultaneously reproduces the second ratio, almost perfectly. 

The huge suppression of strange quark production in PYTHIA would be originated in the absence of correct 
inclusion of strange production in current tunes of PYTHIA. It has been concluded that there is a large increase 
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in the measured production cross section of strange particles as the energy increases from 0.9 to 7 TeV [28j . 
Also, it is found that he difference between predictions of strange particles production and measurements gets 
bigger as the particle mass and strangeness number increase. 
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